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Abstract 

There are presently two models for quantum walks on graphs. The 
"coined" walk uses discrete time steps, and contains, besides the particle 
making the walk, a second quantum system, the coin, that determines the 
direction in which the particle will move. The continuous walks operate 
with continuous time. Here a third model for quatum walks is proposed, 
which is based on an analogy to optical interferometers. It is a discrete- 
time model, and the unitary operator that advances the walk one step 
depnds only on the local structure of the graph on which the walk is 
taking place. This type of walk also allows us to introduce elements, such 
as phase shifters, that have no counterpart in classical random walks. 
Several examples are discussed. 

1 Introduction 

Random walks on graphs are the basis of a number of classical algorithms, 
Examples include, 2-SAT, graph connectivity, and finding satisfying as- 
signments for Boolean formulas. As a result, it is natural to explore the 
quantum counterpart of a random walk, in the hope that it will be useful 
in the development of quantum algorithms. This has led to a number 
of studies. Quantum walks on the line were examined by Nayak and 
Vishwananth 0, and on the cycle by Aharonov, et al. 0. The latter 
study also considered a number of properties of quantum walks on gen- 
eral graphs. Numerical simulations of walks in two and three dimensions 
were performed by Mackay, et al. Absorbing times and probabilities 
of quantum walks on the line were studied by several authors One 
of the main results to come from this work is that quantum walks spread 
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faster than do classical ones. In particular, on the line, the standard de- 
viation of the position of the particle making the walk increases linearly 
with the number of steps rather that with its square root as in the classi- 
cal case. Walks on the hypercube have also been considered, and here the 
results are even more dramatic J! . Kempe has shown that the hitting 
time for the walk from one corner of an n-bit hypercube to the opposite 
corner is polynomial in n for a quantum walk, but exponential for a clas- 
sical one. The quantum walk on the hypercube was subsequently used as 
the basis of a quantum search algorithm The effect of decoherence 
on these walks has also been studied. Brun et al. showed how increas- 
ing dedoherence turns a quantum walk into a classical random walk (3]. 
Kendon and Treganna found that small amounts of decoherence can actu- 
ally speed the convergence of the time-averaged probability distribution 
of a particle in a quantum walk to a uniform distribution |1U) . 

The time steps in the quantum walks considered in these works are 
discrete. Continuous-time quantum walks have also been proposed 
It was shown that on a particular graph, the propagation between two 
properly chosen nodes is exponentially faster in the quantum case. 

There is, at the moment, only one algorithm based on quantum walks 
on graphs, and it was proposed very recently by Childs, et al. |15f. They 
constructed an oracle problem that can be solved by exponentially faster 
on a quantum computer than on a classical one by utilizing continuous 
quantum walks. The vertices of the graph are numbered (named), and 
there are two special vertices called the entrance and the exit. The prob- 
lem is, given the name of the entrance, and the oracle, find the name of 
the exit. The oracle specifies the graph, which belongs to a particular 
set of possible graphs, by taking a binary number as its input, and either 
telling you that this number does not correspond to a vertex, or if it does, 
telling you the names of the adjacent vertices. 

All of the discrete-time quantum walks are based on a particular model, 
the, "coined quantum walk", due to Watrous I16| . In trying to formulate 
a quantum walk on a graph, the most natural thing to do is to let a set 
of orthonormal basis states correspond to the vertices of the graph. If a 
particle is in the state |n), that corresponds to its being located on vertex 
n. Trying to define a unitary evolution using this scheme soon leads 
to serious problems, as was first noted by Meyer (12| . Watrous solved 
this problem by enlarging the Hilbert space in which the quantum walk 
takes place. How this scheme works is most easily seen by considering 
the quantum walk on a line. The vertices are labelled by integers, and, 
in addition, there is a quantum coin, which has two states, \L) and \R), 
corresponding to left and right, respectively. A basis for the Hilbert space 
describing this system is given by the states |n) (g)a), where n is an integer, 
and a is either L or R. A step in this walk consists of applying the 
Hadamard operator, H, to the coin, 

H\R} = -1={\L)-\R)). (1) 
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and then the operator 

Vh = S^\R){R\+S'' (»\L){L\, (2) 
where S is the shift operator, whose action is given by 

S\n) = \n + 1) S''\n) = \n-l). (3) 

The coined quantum walk can be extended in a simple way to regu- 
lar graphs, i.e. those in which all vertices have the same number of edges 
emanating from them. When this is not true, things become more compli- 
cated, and it seems to be necessary to consdier the global structure of the 
graph in defining the walk. So far, no studies of discrete-time quantum 
walks on graphs that are not regular have appeared. 

What we wish to propose here is a different type of discrete-time quan- 
tum walk. It is based on thinking about the graph as an interferometer. 
The vertices are optical elements known as 2A''-ports, where A'^ is the num- 
ber of edges meeting at the vertex, and the edges correspond to paths a 
photon can take through the interferometer. There is no quantum coin in 
these walks. The states are labelled by the edges rather than the vertices 
in the graph, and each edge has two states. If the edge is labelled ab, a 
corresponding to one end and b to the other, then one state is ab, corre- 
sponding to a photon going from a to b, and the other is 6a, corresponding 
to a photon going from b to a. This approach is easily extended to arbi- 
trary graphs; one simply writes down the transition rules for each vertex, 
and all of them taken together define a unitary operator that advances 
the walk one step. In addition, we can add elements to this walk that 
correspond to the addition of phase shifters to paths in an interferometer. 

This model of a quantum walk on a graph is closely related to the 
optical networks considered by Torma and Jex |13|. They considered 
two-dimensional arrays of beam splitters and the propagation of photons 
through them. The horizontal motion of the photon in these networks 
corresponds to the time steps in a quantum walk, and the vertical position 
of the photon is just the position of the particle in the quantum walk. 
Note that these networks provide one with the opportunity to simulate 
the model of quantum walks proposed here with linear optics. 

2 Walk on the cycle 

Perhaps the simplest walk is that on a cycle or ring. Let us label the 
vertices by the numbers through N — 1, where the vertex is identical 
to 0. That is, if we move one step forward from the vertex A'^ — 1, we end 
up at vertex 0. The states of the system are \ j, k), where k = j ±1, which 
can be thought of as a photon on the edge between vertices j and k going 
from j to k. Because each edge has two states, and there are A'^ edges, 
the dimension of this space is 2A'^. 

The vertices can be thought of as beam splitters. Consider what hap- 
pens when a photon travelling in the horizontal direction hits a vertical 
beam splitter. The photon has a certain amplitude to continue in the 
direction it was going, i.e. to be transmitted, and an amplitude to be 
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change its direction, i.e. to be reflected. The beam sphtter has two input 
modes, the photon can enter from eitlier tlie riglit or the left, and two out- 
put modes, the photon can leave heading either right or left. The beam 
splitter defines a unitary transformation between the input and output 
modes. 

We now need to translate this analogy into transition rules for our 
quantum walk. Suppose we are in the state \j — If the photon is 

transmitted it will be in the state \j,j + 1), and if reflected in the state 
\j, J — 1) ■ Let the transmission amplitude be t, and the reflection amplitude 
be r. We then have the transition rule 

\j-iJ)^t\j,j + l)+r\j,j-l), (4) 

where unitarity implies that \t\'^ + |rp = 1. The other possibility is that 
the photon is incident on vertex j from the right, that is it is in the state 
+ If it is transmitted it is in the state \ — and if it is reflected, 
it is in the state \j,j + 1). Unitarity of the beam splitter transformation 
then gives us that 

\j + -i-J)^t*\j,j-l)-r*\j,j + l). (5) 

These rules specify our walk. 

The case t = 1 and r = corresponds to free particle propagation; a 
"photon" in the state \ j, j + l) simply moves one step to the right with each 
time step in the walk. If r 0, then there is some amplitude to move both 
to the right and to the left. A physical system to which this is analogous 
is the motion of a particle in a periodic potential. The beam splitters can 
be thought of as scattering centers with the scattering resulting from a 
localized potential. As is well known, this leads to energy bands, and, as 
we shall soon see, a similar structure emerges in quantum walks on the 
cycle. 

One way of approaching the study of the dynamics generated by this 
walk is to find the eigenvalues and eigenstates of the unitary transforma- 
tion, U, that moves the system a single step. In order to do this, we first 
note that U commutes with the translation operator, T, where 

T\j,j + l) = \j + l,j + 2) T\j + l,j} = \j + 2,j + l). (6) 

This implies that these operators can be simultaneously diagonalized. The 
eigenvalues of T are e'**^ , where 9k = 2iTk/N, and A; = 0, 1, . . . iV- 1. Each 

of these eigenvalues is doubly degenerate, and the two dimensional space 
of eigenvectors corresponding to e**'= is spanned by 

JV-l 

^ 3=0 
N-1 

The eigenstates of U are just linear combinations of \uk+) and \uk-). 
Defining 

IV'fc} = ak+\uk+} + ak-\uk-}, (8) 
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we find that the equation U\tpk} ~ ^li^k) becomes 

T ..-;.■.)( ::!)-(:::)■ <») 

Expressing t as t — |f|e"', we find that the eigenvalues are 

\k± = \t\ cos(6lfc - r?) ± i[l - \t\^ cos'(6lfc - r;)]'/^ (10) 
and the corresponding eigenfunctions are given by 

„(+) ^ !l 

"+ [2Ck{Ck + Sk)V/^ 

(+) ^ -i{Sk + Ck) .^^s 



for Xk+, and 



[2Cfe(Cfe-Sfc)]i/2 



for Afc_. Here we have defined 

= [l-ltl2cos2(efc-r,)]i/2 
Sfc = |t|sin(6»fc -r?). (13) 

One thing we notice immediately, is that for all of these eigenstates, 
the probability to be located on an edge, is the same for all edges, just 
1/A^. That means, that for any initial state, l^'in), the average probability 
distribution, 

m — 1 

= + MU'\^^n}f + \{J + 1, j|t/^|*„oi', (14) 

k = 

where Pj™' is the average probability of being on the edge between j and 
j + 1 after m steps, goes to a constant as m — > oo, if all of the eigenvalues 
in Eq. HUH are distinct [5]. This will be the case if {Nri)/n is not an 
integer. 



3 Walk on the line 

The quantum walk on the infinite line can be approached directly or as 
the limit of the walk on the cycle as N goes to infinity. Because we have 
just found the eigenstates and eigenvalues for the walk on the cycle in 
the previous section, we shall adopt the latter course here. In particular, 
we want to examine what happens when we start the walk in the state 
|0, 1). We shall present numerical results and then follow the approach 
developed by Nayak and Vishwanath to study the long time limit of the 
probability of being on the edge between the vertices j and j + 1. 
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Figure 1: Probability distribution for quantum walk after 50 steps 
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Figure 2: Probability distribution for quantum walk after 1000 steps 

The probability distribution for the particle after n steps can be com- 
puted in a straightforward manner. Wc display the results for the case 
t = r — I/V2 and the initial state |0, 1). In Figure 1 wc have n = 50 and 
in Figure 2, n = 1000. Note that, as with the coined quantum walk, these 
probability distributions are not normal distributions. In addition, the re- 
gion in which the probability of finding the particle is large is, roughly,in 
the case of n = 50, between —35 and 35, and in the case n = 1000, 
between —700 and 700. In both cases this corresponds to the high prob- 
ability region lying between —\t\n and \t\n. This feature of the dynamics 
will be confirmed by our asymptotic analysis. 

The asymptotic probabilities can be calculated as follows. Denoting 
the eigenstates corresponding to the eigenvalues A^-i- and Afc_ as 1^"*;+) 
and lipk-), respectively, we find that the wave function of the particle 
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executing the walk, l^'('r)) is, after r steps, 

JV-1 



k=0 s=± 
JV-1 



^ ^(4+'*Al+|v.+) + 4;)*Ar_"iv.-)), (15) 



where, as mentioned in the previous paragraph, the initial state is taken 
to be |0, 1). The amplitudes to be in the states \ j.j + 1) and \ j + at 
time T are given by 

JV-1 

(i,i + l|*(r)) = ^5]e'^*MAl+|a<+f + Al_l4;f ) 

k = 
JV-1 

(j + l,il*(r)) = ^^e'^'MAl+aif*aW + Al_ai;'-at').(16) 

k = 

The probability of being on the edge between vertices j and j + 1 at time 
T, p{j,j + l;r), is 

pU,J + l;r) = \U,j + l|*(r))|' + \{j + 1, j|*(r))p. (17) 

In order to go the the N ^ oo limit, we need to define a number 
of functions of the continuous variable 9 rather than expressing them as 
functions of the discrete variable Ok- We first define 

C{e) = [l-|tl'cos2(6'-77)]i/^ 

s{e) ^\t\sm{e -7]). (18) 

The eigenvalues also become functions of 0, and we shall express them as 

A±(e) =e'"'±'''', (19) 

where < uj+{0) < n and 

-+(^)=tan-(Mi5|pl), (20) 
and — TT < u!- (9) < and 

We can now proceed to take the N ^ <x limit of the sums appearing in 
Eqs. 1161 . For the sums appearing in the first of these equations we have 



iV — 1 

— 2^6 Ak+\a^^_^\ 



2ir 

^g)g»[j« + ^"+(«)l. 



iV^ Jo 4nC{0)[Ci9) + Sm 

fe=0 

(22) 



'lnC{9)[C{9) ~ S{9)]' 
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The sums in the second equation become 



JV-l 

1 \^ f.^l^kVr „( + )*„( + ) 

k=0 
N-1 

— 2^ e ^k-ak+ flfc- 

fc=0 

We are now going to analyze these integrals in the large r limit by 
using the method of stationary phase. This will be done in two different 
ways. In the first, j will be fixed and r will go to infinity. In the second, 
we shall set j = ar, and then let r go to infinity. Some of the details 
of this analysis are given in the appendix. Here we shall just present the 
results. In the case of fixed j we have that 

+ ~ ^{[l + (-iy+ncos'(r/x + 7r/4) 

■KT\t\ 

+ [l-(-ir+nsin^(T/x + 7r/4)}, (24) 

where < yu < 7r/2, and 

/^ = tan-i(|^)- (25) 

We note that this implies that for any interval located symmetrically about 
the origin, the probability of being in that interval goes like 1/r, whereas 
for a classical random walk starting at the origin, it would go like Xj^/r. 
This implies that, as with the coined quantum walk, this quantum walk 
spreads faster than a classical one. In the case that j = qt, we find that 
there are stationary phase points only if a < \t\. That means that for 
a > p(j, j + 1; r) = p(aT, ar + 1; r) decreases faster than any inverse 
power of r. For a < \t\ we have that p{j,j + 1; r) goes like 1/r. Therefore, 
it is most probable that the particle is located in the region | j| < |i|r, and 
we can say that the allowed region for the particle expands with speed \t\. 



4 Relation between quantum walks 

We now have two different quantum walks on the line, the coined walk, 
where one moves between vertices, and, what we shall call the edge walk, 
where the quantum particle making the walk resides on the edges between 
the vertices. It would be useful to know if the two different walks are 
related. In this section we shall show that they are unitarily equivalent. 
It should be emphasized that this result will only be demonstrated for the 
line, whether it holds for more general graphs is not known. Presently, no 
description of a coined walk for a general graph has appeared. 

Let us begin by examining the Hilbert spaces for the two different 
quantum walks. The canonical ortlionormal basis states of the Hilbert 
space for the coined walk on the line are given by \j £ Z}, 

where the state \j) corresponds to the jth vertex, and \r) and \L) are the 
coin states. The Hilbert space in which this walk takes place is just 
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L'^{Z) ®L^{Z2)- The canonical orthonormal basis of the Hilbert space for 
the edge walk is + 1), |j + 1, G Z}, and the Hilbert space itself 
is L'^{Z X Z2), which is identical to L^{Z) ® L^{Z2). 

Let us now move to the dynamics. The unitary operator, V , that 
advances the coined walk one step is given by 

V = (S® \R){R\ + ® \L){L\){I ® G), (26) 

where G € U{2) is a generalized "coin-flip" operator, and is given by 

G\R} = t\R} + r\L} 

G\L) = -r*\R)+f\L). (27) 

The unitary operator, U, that advances the edge walk one step was given 
in Section 2, and is 

U\j-l,j} = t\j,j + l}+r\j,j-l} 

U\j + l,j) = t*\j,j-l)-r'\j,j + l). (28) 

Define the unitary operator E, which takes L^{Z x Z2) into itself, and is 
given explicitly by 

E\j-l,j} = 

+ = \j)<g>\L). (29) 

We find that 

VE = EU, (30) 

so that at the level of amplitudes, the two walks arc unitarily equivalent. 

There is, however, a difference in the probabilities. In the coined walk, 
the probability to be on vertex j is given by combining (taking the squares 
of the magnitudes and adding) the amplitudes for the states \j) ® \R) and 
\j) Cg> \L). Under the mapping E~^ , these states correspond to states 
on different edges, |j — and \j + respectively. However, the 

probabilities in the edge walk are computed by combining the amplitudes 
for being on the same edge, e.g. those for |j — 1, j) and \ j, Therefore, 
there will be a difference in the probabilitcs for the two walks. This can 
be seen explcitly if we examine the probability distribution for the case 
t = r = l/\/2. We again start in the state |1,0), let the walk go for 50 
steps, but now compute the probability that the particle is on a vertex, 
instead of computing the probability that it is on an edge. The result is 
shown in Figure 3. By comparing this figure to Fig. 1, wo see that the 
overall shape of the probability distributions is similar, but the details are 
different. 

5 Phase shifters 

Going back to the interferometer analogy, we note that wo can add a 
new element to quantum walks that has no analogue in classical random 
walks. Interferometers are made up of multiports and phase shifters; a 
phase shifter imparts a constant phase to a photon that passes through 
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Figure 3: Probability distribution for n = 50, and particle on vertices instead 
of edges. 



it. Suppose we were to put a phase shifter that imparts a phase shift of (j) 
just before the j^^ vertex. The transition rules for the states adjacent to 
this vertex are modified, while the rules for all other states are unaffected. 
In particular, we now have 

+ ^ -r*\j,j + l}+t*e'^\j,j-l}. (31) 

Insertion of a phase shifter into an edge can change the properties of a 
quantum walk, because it changes how difi'crcnt paths interfere. 

One system that allows us to see their effect on the average probability 
distribution is a modified walk on a cycle. Suppose that the number 
of vertices is oven, and that wc put a phase shifter in all of the edges 
whose left end is an even numbered vertex, i.e. every second edge has a 
phase shifter in it. This system is exactly solvable, and by examining its 
eigenstates, we shall see how the average probability distributions it gives 
rise to depend on the value of (j). 

The unitary operator that advances this walk one step acts in the 
following way if j is even 

U2\j,j + l) = te^^\j + l,j + 2)+re'''^\j + l,j} 

U2\j + l,j) = t*\j,j-l)-r'\j,j + l), (32) 

and if j is odd, then 

U2\j,j + l) = t\j + l,j + 2)+r\j + l,j) 

U2\j + l,j} = t*e'^\j,j-l}-r*\j,j + l). (33) 

This operator commutes with translations by two steps, i.e. with the op- 
erator T^. The eigenstates of are given by 

j=0,even 
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j=0,odd 
JV-1 

J— t), even 
/ JV-1 



= ^ e^'^'lJ + hj)- (34) 

j=0,odd 

Each of these states has the eigonvahic exp (—2i6k)- Eigenstatos of U2 are 
just hnear combinations of the above states. In particular, expressing the 
eigenstate of U2, IV*!^^); as 

IVf > = a^Mi) + ak-\ui'l) + bk+\u<-:l) + bk-\ui"l)^ (35) 
the eigenvalue equation U2\ipl^^) = X\tpl^^) becomes 



/ -r* te"'* 

t*e'^'' r 

re^*"* t*e*(*'=+'^) 

\ te'^'^-»k) _r* 



] 
j 



( ak+ \ 




( ak+ \ 


bk- 




bk- 


a-k- 




ak- 


\ bk+ J 




\ bk+ J 



(36) 

The eigenvalues satisfy the equation 

A* + A^[|r|^(l + e^^'*) - e*'*(r^e'''*'= + fe-^'"")] + e^"*" = 0. (37) 

The eigenstates of this system no longer give rise to constant proba- 
bility distributions; the probabilities of being on an even edge (an edge 
whose left-most vertex is even) and an odd edge (an edge whose left-most 
vertex is odd) axe, in general, different. If = 0, these probabilities are 
the same, but if = 7r/2, then this is no longer the case. In the latter 
case we find that 

= i\tf cos(26ifc - 27?) ± [1 - 1*1" cos^(26ifc - 2??)]^/^ (38) 

Choosing the plus sign in the above equation, we find that for an cigen- 
function corresponding io 9k — rj = 7r/4, we have for the ratio of the 
probability of being on an even edge to being on an odd one 

Peven _ 1 + \f\ /on-. 

If Sfc — ?7 cannot be exactly equal to 7r/4 because of the values of N or 
ri, then for 9k — v close to n/4 the ratio of even-edge to odd-edge prob- 
abilities will be approximately given by the above equation. This ratio 
is not generally equal to one, which means that the average probability 
distribution to which a given initial state converges will not be constant. 
The introduction of the phase shifters has changed the character of the 
quantum walk. 

These changes can also be seen by calculating the probability distri- 
butions after a finite number of steps. This is done for 50 steps and for 
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Figure 4: Probability distribution for n = 50 and (p = Tr/2. 



the case t = r = l/\/2 and (f) = 7r/2, 7r/3 in the following figures. The 
initial state is, as before, |0, 1). These can be compared to Fig. 1, which 
corresponds to the case (p — 0. It can be seen that the introduction of 
the phase shifter greatly changes the character of the probability distri- 
bution. Note that particularly for the case of = 7r/2, the size of the 
region in which it is very likely that the particle will be found is smaller 
than when </> = 0. For a small number of steps, it is easy to verify by hand 
that destructive interference in the 4> = n/2 case makes the walk spread 
more slowly than when (p = 0, and the numerical results indicate that this 
feature persists for at least 50 steps. 

If we extend this walk to the infinite line, the difference caused by the 
phase shifters can be seen in the asymptotic behavior. For both <f) — and 
(f> = n/2 the size of the region in which it is most likely to find the particle 
grows linearly with the number of steps, but the "speed" is different. We 
saw that in the case = the probability distribution p(j.,j + 1 : r) falls 
off rapidly for \j\ > \t\T. If = tt/2, it falls off rapidly for \j\ > |tpr, 
which means that the size of the high probability region is smaller in this 
case. We see yet again, that phase eff'ects, which do not exist in classical 
random walks, can significantly influence the behavior of quantum walks. 



6 Probability current 



In standard quantum mechanics, it is possible to define a probability cur- 
rent density. In one dimension, if the wave function of the particle is 
'tp{x,t), then the probability current density is given by 



3{x,t) = 



2im 



^*ix,t)^ip{x,t) -■tp{x,t)-^tp*{x,t) 



(40) 
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Figure 5: Probability distribution for n = 50 and (f> = 7r/3. 

where m is the mass of the particle, and we are using units in which h = 1. 
This current has the property that 

^ / ' m^,t)f = -~[j{x2,t)-j{x^,t)], (41) 

that is, the change in the probability of the particle being in a particular 
region is given by the net flow of probability into the region. For an eigen- 
state of the Hamiltonian, the probability density \ip{x,t)\'^ is independent 
of time, so that the probability current density is a constant. We would 
hke to show that there is a quantity similar to the probability current 
density for quantum walks. 

Suppose that the state of the walk on the cycle is given by 

N-l 

1*) = j + 1) + c, + i,,|j + 1, j)). (42) 

3 = 

Define the probability current at the k^^ vertex to be 

..Mc:„.c:_.,,(i.f. _;.)(::;::). m 

We find that, if APk^k+i is the change in the probability of being on the 
edge between vertices k and fc + 1 in one step of the walk, then 

APk.k+i = Jfe+i - Jfc, (44) 

which is the discrete analogue of Eq. 14111 . For an eigenstate of the walk, 
this current, Jk, must be independent of k. The above equation holds even 
if the transmission and refiection amplitudes are different for each beam 
splitter, i.e. they depend on k, and also if phase shifters are present. In 
that case, if tk is the transmission amplitude at vertex k, rt the reflection 
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amplitude, and 4>k the phase shift of the phase shifter just to the left of 
vertex k, the current at this vertex is give by 

Jk = (4+1.. 4-1..) ( I'J^ ) ( 'cu-i ) ■ ^^'^ 

We can use this fact to demonstrate a general property of eigenstates 
of certain kinds of walks on a line. Suppose that all of the beam splitters 
located at vertices fe < and k > N have transmission amplitude t = 1 
and reflection amplitude r — 0. The beam splitters for < k < N can 
have any value of transmission amplitude and reflection amplitude, and 
these values can vary from vertex to vertex. We shall refer to the vertices 
between and N as the scattering region. The problem we are considering 
is analogous to the scattering of a particle moving in one dimension off of 
a potential, which is nonzero only in a bounded interval. From the point 
of view of a walk, we may be interested in a walk that starts to the left 
of the scattering region in a right-moving state, and finding out how long 
it takes to get through the scattering region. 

The eigenstates of this type of walk are of two types. The first of a 
particle coming in from the left, its reflected amplitude, and a transmitted 
amplitude to the right of the scattering region. The second consists of a 
particle coming in from the right, its reflected amplitude, and a transmit- 
ted amplitude to the left of the scattering region. We shall consider the 
first type, which can be expressed as 

JV-l oo 

j = — CXD j = N 

Setting the eigenvalue. A, equal to exp(— iS), the equation Ul"^) = \\'^) 
gives us that 

Ci.^+1 = e'(^+i>'c_i.o forj < -1 

Cj+i.j = e"'^^"''^^''co.-i forj < -1 

Ci,j+i = e''^-'^"cjv,]v+i forj > A^. (47) 

The amplitude c-i,o can be thought of as the amplitude of the incoming 
wave, co,-i the amplitude of the reflected wave, and cn,n+i the amplitude 

of the transmitted wave. We can find a condition that these quantities 
must satisfy, if we make use of the fact that J_i = Jn+i- This gives us 
that 

+|cjv,Af+i| = |c_i,()| , (48) 

where we have used the fact that |c_2.-ip = \c-i,u\^ ■ Defining the re- 
flection coefficient of the scattering region to be 7i = |co,-ip/|c_i.oP and 
the transmission coefficient to be T = |cjv.Ar+i |^/|c-i.o|^, then the above 
equation can be expressed as i? -|- T = 1. 



7 Vertices with more than two edges 

So far we have only considered vertices at which two edges meet, but if 
we are to construct graphs more complicated than lines, we need to see 
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how a vertex with more that two edges emanating from it behaves. We 
shall look at two examples, one with three edges, and another with an 
arbitrary number. 

A vertex with three edges emanating from it, which is inspired by the 
optical multiport known as the fritter I14| . can be described as follows. 
Let us label the vertex with three edges meeting at it by O, and the 
opposite ends of the edges by A, B, and C. The ingoing states for this 
vertex are |^0), \BO), and \CO), and the outgoing states are |OA), \OB), 
and \OC). Setting z = exp(27ri/3), we have for the transition rules 

\AO) 
\BO) 
\CO) 

This vertex has the property that an incoming particle is equally likely to 
exit through each edge. However, note that because the incoming states 
from different edges behave differently with regard to their phases, the 
use of this vertex requires the labelling of edges. In this particular case, 
only one of the edges needs to be labelled. If we attach a label to AO, 
we interpret it to mean that if the input state is along either of the other 
two edges, then the output with a phase factor of z* is along the labelled 
edge. For any edge, if the input state is along this edge, the part of the 
output state along the same edge has a phase factor of one. 

It is possible to define vertices that do not require the labelling of 
edges, though the cost is that the probabilities of exiting through each of 
the edges are no longer the same. This particular kind of vertex is very 
closely related to the quantum coin used in the walk on the hypercube 
[S] |7). Let the vertex at which all of the edges meet be labelled by O, 
and the opposite ends of the edges be labelled by the numbers 1 through 
n. For any input state, \kO), where k is an integer between 1 and n, the 
transition rule is that the amplitude to go the output state \Ok) is r, and 
the amplitude to go to any other output state is t. That is, the amplitude 
to be reflected is r, and the amplitude to be transmitted through any of 
the other edges is t. Unitarity places two conditions on these amplitudes 

{n - + \r\^ = 1 

{n-2)\t\^ + T*t + t*r = Q. (50) 

As an example, for the case n = 3, possible values of r and t are r — —1/3 
and t — 2/3. Because each of the edges in this vertex behaves in the same 
way, they are equivalent to each other and no labelling is necessary. 

In order to construct a walk for a general graph, one chooses a unitary 
operator for each vertex, i.e. one that maps the states coming into a vertex 
to states leaving the same vertex. One step of the walk consists of the 
combined effect of all of these operations; the overall unitary operator,?/, 
that advances the walk one step is constructed from the local operators 
for each vertex. Explicitly, the edge state \ab), which can be interpreted 
as going from vertex a to vertex b, will go to the state Ub\ab) after one 



-^{\OA) + \OB) + \OC)) 

-^{z*\OA} + \OB) + z\OC}) 
v3 

-^{z''\OA)+z\OB) + \OC)). 



(49) 
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step, where Ub is the operator corresponding to vertex b. This prescription 
guarantees that the overall operation is unitary, in particular, U acting on 
any other edge state \cd) will give a state orthogonal to U\ab). If d = 6, 
then \ab) and \cd) will be mapped onto the same set of states (the states 
leaving vertex b), but the unitarity of Ui, will ensure that U\ab) and U\cd} 
are orthogonal. If d 7^ 6, then U maps \ab) and \cd) onto different sets 
of states, and the so results are then orthogonal. Therefore, as the edge 
states make up an orthonormal basis of the Hilbert space in which the 
walk occurs, and U maps this basis to another orthonormal basis, it is 
unitary. 

8 Conclusion 

Quantum walks on graphs seem promising for the development of algo- 
rithms, because they spread over a graph faster than does a classical walk, 
and can thereby explore the structure of the graph faster than can a clas- 
sical random walk. Here we have discussed a discrete quantum walk that 
is based on an analogy to optical interferomenters. The vertices act as 
optical multiports and phase shifters can be inserted into the edges. As 
we have seen, the behavior of this type of walk depends on both types 
of elements. The advantage of this type of quantum walk is that it can 
easily be defined for any graph. 

In this paper, we have, for the most part, confined our attention to 
the quantum walk on the line. We found that the probability distribu- 
tion of the particle making the walk spreads linearly with the number 
of steps, and with a "speed" given by In addition, there is a prob- 
ability current, whose "divergence" gives the probability flowing into an 
edge. This allowed us to define refiection and transmission coefficients for 
one-dimensional graphs. 

The extension of these results to more general graphs is clearly the 
next step. It has been shown how definine a quantum walk for any graph, 
but the properties of these walks have yet to be explored. 
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Appendix 

Here we want to explicitly show how the asymptotic analysis on the in- 
tegrals in section HI. We shall show how to do the analysis for the first 
integral in Eq. 12211 . and the procedure for the others is similar. 

We shall first consider the case when j is fixed. Then the stationary 
phase points are the solutions of the equation = 0, and we find the 

two solutions 6 = rj and 9 = rj + tt. We also find that ioKri) = \t\/\r\ and 
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w"(?7 + it) — —\t\/\r\. In addition, w+{'q) = ^ and to+irj + vr) = tt — /x. 
Inserting these values into the standard formula for stationary phase |17|. 
we find that 

^g,gi[je+r^+(9)] \r\ ^ 1 ( 1^ 



47rC(9)[C(6l) + 5'(6i)] 2 y27rr|t| 

Let us now consider the case when j — ar. The stationary phase 
points are now given by the solutions of 

/ (a-, jtj sin(6' - J]) 

"+(^)^ [l-|t|^cos^(g-,)]V2 =-"- (52) 

For this equation to have any solutions, we find that a < \t\. If this 
condition is not satisfied, there are no stationary phase points, and the 
integral decreases faster than any inverse power of r. If a < then there 
are two solutions, and they satisfy the conditions n < 9 — r/ < 2tt and 

Explicitly, if 7 lies between and n/2 and satisfies 

then the two solutions to Eq. 15211 arc 61 =77 + 7 + 71 and O2 — ri + 2n — "f. 
We find that 

\r\ 

Lj'iie^) = ±{\t\' -a'f/\l-a^), (55) 
and uj+{9i) = tt — 1/ and ^+(62) — v, where < < 7r/2 and 

Finally for the integral we find that 



2ir 



'iTxC{e)[c{e) + s{e)] 



2(1 -a) 



27rr(|tl2 -a2)i/2 



_^^iT[a(2-K — /) + i/]+i7r/4j ^^rj^ 



Finding the asymptotic form of the other integrals in a similar fashion, 
we have for j — ar and a < \t\ 

\r\ 

P{jj + i;r) ~ 77-7^ 5u727i r[l + a(-l)^cos(7rar)] 

{l + (-l)^sin[2r(a7-!y)-7rar]}. (58) 
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